Ij  '  A0»A114  720  LEHIGH  UNXV  BETHLEHEM  PA  CENTER  FOR  THE  APPLICATION— ETC  F/G  20/11 

I  BIFURCATION  IN  A  ELASTIC  PLATE  ON  A  RIGID  SUBSTRATE. <U) 

I  JAN  82  R  S  RIVLIN  N00014-76-C-0235 


UNCLASSIFIED  CAM-100-37  NL 

I 


;  FILE  COPY 


C  ; 

'  i 
r-  ‘ 
r  < 


Lehigh 

University 


BIFURCATION  IN  AN  ELASTIC  PLATE  ON  A  RIGID  SUBSTRATE 

by 

R.S.  Rivlin 


TECHNICAL  REPORT  NO.  CAM-100-37  January,  1982 

OFFICE  OF  NAVAL  RESEARCH  CONTRACT  NO.  N00014-76-C-0235 


82 


Center  for  the 
Application  of 
Mathematics 
(CAM) 

05  21  03  6: 


Bifurcation  in  an  Elastic  Plate  on  a  Rigid  Substrate 

by 

R.S.  Rivlin 

Lehigh  University,  Bethlehem,  Pa. 


Abstract 

An  infinite  plate  of  neo-Hookean  elastic  material  is 
bonded  on  one  face  to  a  rigid  substrate.  It  is  subjected  to 
a  uniform  shear  and  dead- loaded  with  a  uniform  thrust.  A 
periodic  bifurcation  solution  is  obtained  when  the  thrust  per 
unit  area  exceeds  a  critical  value.  The  relation  between  the 


wave-length,  thrust,  amount  of  shear  and  plate  thickness  is 


obtained. 
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1.  Introduction 

In  this  paper  we  consider  the  critical  loading  conditions 
for  which  a  bifurcation  solution  is  obtained  when  an  infinite 
plate  of  incompressible  isotropic  neo-Hookean  elastic  material, 
bonded  on  one  face  to  a  rigid  substrate,  is  subjected  to  a 
uniform  shear  of  amount  K  and  simultaneously  dead- loaded  on 
the  other  face  by  a  uniform  normal  thrust. 

We  suppose  that  an  infinitesimal  spatially  periodic  defor¬ 
mation  is  superposed.  A  secular  equation  is  obtained  for  the 
determination  of  the  wave-length  of  the  superposed  deformation. 
This  secular  equation  yields  a  real  value  for  the  wave-length 
when  some  critical  value  of  the  normal  thrust  per  unit  area, 
which  depends  slightly  on  K  ,  and  is  approximately  equal  to 
twice  the  shear  modulus,  is  reached.  For  this  critical  value 
of  the  thrust,  the  wave-length  of  the  infinitesimal  superposed 
deformation  is  zero  and,  as  the  thrust  is  increased  beyond  this 
value,  K  remaining  fixed,  the  wave-length  increases. 

Except  at  values  of  the  thrust  per  unit  a_*ea  near  the 
critical  value,  the  wave-length  is  proportional  to  the  square 
root  of  the  thrust  per  unit  area  and  is  nearly  independent  of 
K  .  At  values  of  the  thrust  near  the  critical  value,  the  wave¬ 
length  becqftea  e>  ly  sensitive  to  the  thrust.  At  all 
values  of , the^  thu’us  .  i  of  K  ,  the  wave-length  is  propor¬ 
tional  to  tfte.  thickness  of  the  plate,  as  may  be  expected  from 

dimensional  considerations. 

,  • 

For  specified  values  of  the  thickness  and  of  the  thrust 

per  unit  area,  beyond  the  critical  value  of  the  latter,  a 

«  • 

periodic  static  bifurcation  solution  is  obtained  with  uniquely 


determined  wave-length.  However,  if  the  plate  were  finite  in 
the  direction  of  the  periodicity,  the  end  conditions  would  en¬ 
able  us  to  determine  the  spectrum  of  values  of  the  wave-length, 
and  hence  of  the  thrust  per  unit  area,  at  which  the  bifurcation 
solution  can  occur. 


2. 


Basic  equations 

We  consider  a  flat  plate  of  incompressible  isotropic  neo- 
Hookean  elastic  material  of  thickness  h  to  be  located  with 
its  major  surfaces  normal  to  the  2 -axis  of  a  rectangular  Cartesian 
coordinate  system  x  .  The  dimensions  of  the  plate  in  the  1  and 
3  directions  are  supposed  large  compared  with  h  . 

Let  5  be  the  vector  position  of  a  generic  particle  of  the 
plate  in  its  undeformed  state  and  let  £a(a»l,2,3)  be  the  com¬ 
ponents  of  €  in  the  system  x  .  We  suppose  that  initially  the 
plate  is  located  with  its  major  surfaces  in  the  planes  C2  *  0 
and  C2  *  h  and  that  the  face  52  -  0  is  bonded  to  a  rigid 
plate  which  remains  fixed  in  space. 

Suppose  the  elastic  plate  undergoes  a  deformation  in  which 
a  particle  initially  at  £  moves  to  vector  position  x  with 
components  x.(i«l,2,3)  in  the  system  x  .  Then,  the  strain- 
energy  W  per  unit  volume  is  given,  in  appropriately  chosen 
units,  by 


W 


7Cxi 


x . 


i,a  i,a 


.-3) 


(2.1) 


where  ,a  is  used  to  denote  differentiation  with  respect  to  £a  . 
The  Piola-Kirchhoff  stress  ff  j  ,  referred  to  the  system  x  , 
is  given  by 


ai 


i*ot 


7  PeijkEa$YXj,6xk,y  » 


(2.2) 


where  e.^  denotes  the  alternating  symbol  and  P  is  an  arbi 
trary  hydrostatic  pressure.  Since  the  material  considered  is 


5. 


incompressible 


det  lx 


1.0 


(2.3) 


We  now  suppose  that  the  deformation  £  -*■  x  is  the  resul- 

*v  ** 

tant  of  a  finite  deformation  K  X  and  an  infinitesimal  defor 

<*»  «w 

mat  ion  X  -*•  x  ,  where 


x  -  X  +  eu  (2.4) 

*» 

and  e  is  a  small  parameter.  We  assume  that  the  force  system 
associated  with  the  deformation  £  -*■  x  differs  by  terms  of 

order  e  from  that  associated  with  the  deformation  £  -►  X  . 

*«•  •*» 

We  accordingly  write 


"ai  "  nai  +  eiToi  »  p  ■  p  +  eP  •  (2.5) 

Then,  from  (2.2)  we  obtain,  with  (2.4)  and  (2.5), 

nai  "  Xi,a  '  1  PEijkEa0yXj . SXk, y  * 

woi  "  Ui»o  -  I  Eijk£aey{P(Xk,yuj,B+  Xj.8uk,y> 

+  pXj.SXk,y>  ‘  <2*6> 


Similarly  from  (2.3)  we  have 


detlXi,«l  "  1  »  eijkEagrxi,aXj,Buk,y  "  0  *  (2,7) 


6. 


The  Piola-Kirchhoff  equations  of  equilibrium  yield 


n 


ai,a 


0  and 


ffoi,a 


0  . 


(2.8) 


7. 


3.  The  governing  equations 

If  the  deformation  5  +  X  is  a  simple  shear  of  amount  K  , 

<v  <v  * 

for  which  the  direction  of  shear  is  the  1-direction  and  the 
plane  of  shear  is  the  12 -plane,  then 


X1  "  Ci+  U2  *  X2  *  C2  »  X3  ‘  C3  *  (3.1) 

We  shall  assume  that  the  superposed  infinitesimal  deformation 
X  -*■  x  is  a  plane  deformation  in  the  12-plane.  We  can  then  write 


U1  *  Ul^l»C2)  *  U2  “  U2C51,52J  »  u3  "  0  .  (3.2) 

With  (3.1)  equation  (2.7)1  is  automatically  satisfied  and 
with  (3.1)  and  (3.2)  equation  (2.7)2  yields 


ul,l  +  u2,2 


Ku 


2,1 


(3.3) 


Also,  with  (3.1),  (3.2)  and  (3.3)  equations  (2.6)  yield 


nn  •  n22  -  n33  -  i  -  p  ,  n12  »  kp  ,  n21  -  k  , 


n31  "  ni3  "  n23  n32  "  °  * 


(3.4) 


and 


*11  •  Ul,l  •  Pu2 , 2  *  P  •  *22  '  u2,2'  Pul,l-  P>  *3,3’"P> 

(3.5) 

*12*U2,l*Pul,2*plt>  *21-ul,2*Pu2,l-*3l’*13"*23"*32*  °- 


8. 


With  (3.4),  equation  (2.8)^  implies  that  P  is  constant. 
With  (3.5),  the  incremental  equation  of  equilibrium  (2.8)2 
yields 


u  +  u 

1,11  1,22 


u 


2 , 11+  U2 ,22 


+  »1“  Pi 


*3 


0.  (3.6) 


The  last  of  these  equations  implies  that  p  ■  p(C1,52)  . 

From  (3.3)  it  follows  that  there  exists  a  function 
<K51»C2)  in  terms  of  which  we  may  express  1^,  u2  by  the 
relations 


U1  "  ^ * 2  "  K*»l  *  U2  "  “  *»!  • 


(3.7) 


By  substituting  (3.7)  in  (3.6)2  2  ®nd  eliminating  p  ,  we  obtain 


V2[(1+K2)*,i;l  +  *,22  -  21*.  12]  -  0  ,  (3.8) 


where 


V 


2 


(3.9) 


We  shall  obtain  solutions  of  (3.8)  which  are  sinusoidal 
in  the  1- direction.  Accordingly,  with  the  usual  complex  nota¬ 
tion,  we  write 

4>  -  ♦U2)e  1  ,  (3.10) 

where  k  is  a  constant.  Introducing  (3.10)  into  (3.8),  we 
obtain 


9. 


f  •  *  «  ftt  o  9  t*  *  »  li  O 

<fr  -  2ikK$  -  k2(2+lT)$  +  2ikJK*  +  k4(l+K2)$  -  0  ,  (3.11) 

where  the  prime  denotes  differentiation  with  respect  to  C2  . 
Equation  (3.11)  has  the  general  solution 


°A«2 


r  c. 

<f>  -  l  a.e  , 


(3.12) 


where  the  a’s  are  (complex)  constants  and  the  a's  are  given 
by 


a  -  k  ,  a  -  -k  ,  a  «  k(l+iK)  ,  a.  -  k(-l+iK)  . 


With  (3.10)  and  (3.12),  we  obtain  from  (3.7) 


(3.13) 


ikCi  k  a.g 

•  1  i  (oA-ikK)aAe  A 

A-l  A  A 

./*«,  . 

A-l  A 


(3.14) 


With  (3.14),  equations  (3.6)  yield 


ik£,  4  2  9 

P»x  ■  e  1  Z  (°A"k2) («A-ikK)aAe 


aA52 


ike,  it  9  9  9  aACj> 

P» 2  "  •  1  I  (o|-k2)[-ik(l+K2)  ♦  KoA]aAe  A  2  . 


(3.15) 


These  equations  yield 


P  -  -  i  elkC\i!_  (o2-k2)(oA-xkK)aAeaAC2 


♦  constant 


(3.16) 


With  (3.14)  and  (3.16),  we  obtain  from  (3.5) 


I 

IT 


V 


IT 


11 


22 


12 


IT 


21 


e 


e 


e 


lkEi  i 

A-l 

A-l 

A-l 

A-l 


^.k«A(p  ♦  £|)  ♦  . 

(-HcaA(2  ♦  ^^-k2-k2P)]aAeV2 

a2 

{-iklCaA(p-l  +  -A) 

+  k2(l+K2)  ♦  a2(P-K2)|  aAe°A^2 

^-ikKaA+(oi2+k2P)jaAeaA52  .  . 


j 

i 


11. 


4.  Solution  of  the  secular  equation 

In  this  section  we  shall  assume  that  the  surface  conditions 
on  the  surface  52  ■  h  are  dead-loading  conditions.  Accordingly, 


w22 1  C2=  h  =  "21 U2-  h  “  0 


(4.1) 


On  the  surface  ■  0  ,  we  have 


UlU2  =  0  =  u2  €  -  0  "  0  ‘ 


(4.2) 


Introducing  (3.14)  into  (4.2),  we  obtain 


l  (aA-ikK)a.  -  0  ,  l 


a. 


A=  1 


A-l 


Again  introducing  (3.17)2  4  *nt0  (4.1),  we  obtain 


k  k 

l.  BAaA  *  0  *  .1-  CAaA  •  0  * 


A-l 


A-l 


where 


B, 


2 

^-lkaA(2+P-  K(a2.k2.k2p)j 

^-tkKaA  +  aA+k2pJeaAh  . 


aAh 


(4.3) 


(4.4) 


(4.5) 


The  necessary  and  sufficient  condition  for  (4.3)  and 
(4.4)  to  have  a  non-trivial  solution  for  aA  is 


12. 


(4.6) 


With  the  expressions  (4.5)  for  and  ,  the  secular 

equation  (4.6)  can  be  rewritten  (see  Appendix)  as 


B^2  +  2B2Q  ♦  B3  *  0  , 

where 

Q  »  1  ♦  P  , 

B  -  (K2+4) - (K2cosh2u  +  4cos  Ky)  , 

$2  *  2K2(cosh2y-cos  Ky)  , 

B3  -  K2[(K2+4)  ♦  (K2cosh2y  +  4cos  Ky)]  , 

with 


(4.7) 


(4.8) 


y  ■  kh  . 

From  (4.7),  Q  is  given  by 


(4.9) 


(4.10) 


We  note  from  (4.8)  that 


13. 


$2-ei63  ■  K2(K2+4)(K2sinh22y-4sin2Ky)  , 

61  -  -2 (K2sinh2y-4sin2yKy)  ,  (4.11) 

&3  -  2K2(K2cosh2v»+4cos2jKy)  . 

Since  K2sinh22y  >  4sin2Ky  for  all  K  and  y  ,  it  follows 
from  (4.10)  and  (4.11),  that  Q  is  real  for  all  K  and  y  . 
Also,  it  is  evident  from  (4.11)2  3  that  B1  <  0  and  >  0  . 

It  follows  that  the  negative  alternative  in  (4.10)  leads  to  a 
negative  value  for  Q  .  Since,  from  (3.4),  the  condition  for 
the  normal  traction  on  C2  ■  h  to  be  a  thrust  is  P  >  1  ,  i.e. 

Q  >  2  ,  it  follows  that  the  negative  alternative  in  (4.10) 
corresponds  to  the  normal  traction  on  -  h  being  tensile. 

We  note  also  that  $2  0  .  Accordingly,  the  necessary 

and  sufficient  condition  that  Q  >  2  ,  i.e.  the  normal  traction 
is  a  thrust,  is 

(^'‘W’5  >  *  (4‘12) 

From  (4.8)  it  follows  that 

2B1+  S2  "  2(K2+4)  (1-cos  Ky)  .  (4.13) 

Accordingly,  apart  from  the  trivial  case  when  cos  Ky  *  1  , 
the  inequality  (4.12)  is  always  satisfied  by  the  positive  alter¬ 
native  in  (4.10).  We  conclude  that  this  corresponds  to  the  nor¬ 
mal  force  on  the  surface  52  ■  h  being  a  thrust.  In  the  remain 
der  of  this  paper  we  consider  only  this  situation. 


In  Fig.l  we  plot  from  (4.10),  with  (4.8)  and  (4.11),  the 
values  of  the  thrust  N  »  P-l-Q-2  against  y  for  K  fixed. 
The  curve  in  Fig.l  was  drawn  for  K  -  0.4.  However,  the  depen 
dence  of  N  on  K  ,  for  K  in  the  range  0  to  0.4,  is  less 
than  2%.  It  was  found  from  the  numbers  from  which  the  figure 
was  plotted  that  Q  -  2  a  13.2/y2  for  y  <  1  . 

It  is  instructive  to  consider  three  limiting  cases. 

(a)  y  fixed.  K  ■»  0  . 

From  (4.8)  it  follows  that 


61  -  2K2(y2-sinh2y)  +  0(KU)  , 

82  -  4K2sinh2y  +  0(KU)  , 

B3  -  8K2  +  0(KU)  . 

With  these  expressions  we  obtain  from  (4.10) 

2{sinh2y  +  (sinh^y  +  sinh2y  -  y2)*5} 

Q  "  o  2  * 

sinh  y  -  y 

In  the  limit  y  •  ,  Q  ■  4  . 

In  the  limit  y  *  0  ,  Q  -  ^-(1+  -)  . 

y  /J 

(h)  K  fixed,  y  +  *  . 


(4.14) 


(4.15) 


From  (4.8),  it  follows  that,  as  y  ♦ 


9 


15. 


Then,  from  (4.10),  we  obtain 
Q  -  2  +  (K2+4)’5  . 

As  K  increases  from  0  to  1  ,  Q|y-—  increases  from 

4  to  4.236. 

(c)  K  fixed,  u  -*•  0  . 

From  (4.8),  we  obtain,  as  y  *►  0 

61  "  -^U|c2Ck2+4)  +  °Cy6)  , 

62  -  y2K2(K2+4)  +  0(yU)  , 

83  -  2K2 (K2+4)  +  0(y2)  . 

Then,  from  (4.10),  we  obtain 

Q  -  ^  d+^y)  *  12.93/y2  , 

in  agreement  with  the  result  obtained  in  case  (a) . 


(4.17) 


(4.18) 


(4.19) 
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5.  Appendix 

With  the  expressions  (3.13)  for  aA  ,  equation  (4.6)  can 
be  written  as 

CBlVB2Cl>  ♦  tB2C3-B3C2)  *  CB3Ct-BliC3)  *  (»hCx-B1CkJ 

♦  ^K{(B1-B2)(C3-Cll)  -  (C1-C2)  (B3-BkJ }  -  0  .  (5.1) 

Also,  with  (3.13)  and  the  notation  Q  -  1  +  P  and  y  -  kh  , 
we  obtain  from  (4.5) 

Bx  -  k2[K(l-Q)-tQJey  , 

B2  -  k2[K(l-Q)+iQ)e_y  , 

B3  -  -k2[iQ+K(l+iK)]eCl+lK)y  , 

B.  -  k2[iQ-K(l-iK)]e'C1'lK)y  , 

(5.2) 

Cx  -  k2(Q-iK)ey  , 

C2  -  k2(Q+iK)e’y  , 

C3  -  k2(Q+iK)eC1+'K)y  , 

Cu  -  k2(Q-iK)e‘Cl'lK)y  . 

From  (5.2)  we  obtain 

B3CU  -  *kC2  -  (B^g-BgC^e21^  -  -2xk4{Q  +  (Q-l)K2}e2lky  , 

B2C3  ’  B3C2  "  k4(Q+iK){2xQ+K(l-Q)  ♦  K(l+xK)}elICy  ,  (5.3) 

BUC1  -  B^  -  k4(Q-iK){2iQ-K(l-Q)-K(l-iK)}elKy  . 


17. 


Whence , 

(B2C3-B3C2)  +  CBUC1-B1CU)  -  4ik4(Q2+K2)elKy  .  (5.4) 

We  also  obtain  from  (5.2) 

(b1-b2)(c3-c4)  -  (C2-C2) (B3-B*) 

-  kU{[QK(4-Q)  +  K3] (e2y+e"2y)  +  2K(Q2+K2) }elKy  .  (5.5) 

With  (5.3),  (5.4)  and  (5.5),  the  secular  equation  (5.1)  may 
be  rewritten  as 

-2{Q2+(Q-l)K2}(l+e2lKy)  +  4(Q2+K2)elKy 

♦  jK2{[Q(4-Q)+K2] (e2y+e‘2y)  +  2 (Q2+K2) }elKy  -  0  .  (5.6) 

From  this  the  relation  (4.7)  with  (4.8)  is  easily  obtained. 


\ 
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